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Abstract 

Approximations based on two-particle irreducible (2PI) effective actions (also 
known as <i>-derivable, Cornwall-Jackiw-Tomboulis or Luttinger-Ward function¬ 
als depending on context) have been widely used in condensed matter and non¬ 
equilibrium quantum/statistical held theory because this formalism gives a ro¬ 
bust, self-consistent, non-perturbative and systematically improvable approach 
which avoids problems with secular time evolution. The strengths of 2PI approx¬ 
imations are often described in terms of a selective resummation of Feynman 
diagrams to infinite order. However, the Feynman diagram series is asymp¬ 
totic and summation is at best a dangerous procedure. Here we show that, at 
least in the context of a toy model where exact results are available, the true 
strength of 2PI approximations derives from their self-consistency rather than 
any resummation. This self-consistency allows truncated 2PI approximations 
to capture the branch points of physical amplitudes where adjustments of cou¬ 
pling constants can trigger an instability of the vacuum. This, in effect, turns 
Dyson’s argument for the failure of perturbation theory on its head. As a result 
we hnd that 2PI approximations perform better than Fade approximation and 
are competitive with Borel-Pade resummation. Finally, we introduce a hybrid 
2PI-Pade method. 
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1. Introduction 

Two-particle irreducible (2PI) effective actions and approximation schemes 
based on them are often touted as useful techniques when it is necessary to 
go beyond the standard perturbative field theory, with applications to thermal 
and non-equilibrium plasmas/fluids, strongly coupled quantum field theories 
and systems dominated by many-body collective effects. The technique was 
originally developed by Lee and Yang [1], Luttinger and Ward [2], Baym [3] 
and others in the context of many-body theory, then extended by Cornwall, 
Jackiw and Tomboulis |3| to relativistic field theory where it found its natural 
formulation in terms of functional integrals. Since then a broad literature has 
developed surrounding 2PI effective actions and their generalizations (see [5] for 
a good introductory review). 

Approximations based on 2PI effective actions are often justified as a se¬ 
lective re-summation of perturbation theory (some recent examples: EllZlIHj; 
interestingly though, this point of view is not found in the Cornwall, Jackiw 
and Tomboulis paper 0), however they are not really a summation method in 
the same sense as, for example, Borel summation. Rather, this interpretation 
comes in a rather roundabout way. First, a set of self-consistent equations of 
motion is derived for the mean field ip = {4>) and connected correlation function 
A = ((/>^) — starting directly from the non-perturbative definition of the 
theory through the path integral. Only after formally solving the equations 
of motion by repeatedly iterating them does one obtain the usual perturbative 
expansion for these quantities or, if the equations of motion are truncated, a 
selective re-summation of perturbation theory appears. 

In this work we attempt to clarify the connection between 2PI approxima¬ 
tions, traditional perturbation theory and re-summation methods with a special 
emphasis on what analytic features of the 2PI formalism allow improved approx¬ 
imations to be obtained in the presence of a divergent Feynman diagram series. 
Unfortunately robust comparisons are difficult because the large order behaviour 
of perturbation theory is known, at best, in a sketchy form for most field theories 
of interest. To that end we restrict attention to a genuinely trivial model “field 
theory” in zero spacetime dimensions (i.e. probability theory) for which exact 
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results are easily obtainable and all complications due to renormalization etc. 
disappear. This model nevertheless accurately represents the typical combina¬ 
toric structure of large order perturbation theory, at least in those cases where 
the behaviour is known in more interesting field theories. We also introduce 
a “spectral function” representation of the Green function (similar to the one 
first introduced by Bender and Wu laiini) to capture the non-analyticity of the 
solutions in the various methods. 

The existence of the spectral representation is connected to the branch cut 
of physical amplitudes on the negative coupling (A) axis. This branch cut is 
due to the non-existence of the theory at negative couplings: the path integral 
diverges due to a potential unbounded from below. In a higher dimensional 
field theory this has a simple physical interpretation: the vacuum is unstable 
and, after tunneling through a barrier, the system rolls down the potential m- 
For weak coupling the semi-classical approximation is valid and the tunneling is 
exponentially suppressed, giving an imaginary contribution ~ exp (—1/A) to the 
vacuum persistence amplitude which is inherited by the Green function. This 
exponential behaviour can be seen in the spectral function we obtain. 

Dyson |12j argued that a very similar phenomenon occurs in quantum elec¬ 
trodynamics (QED). We briefly reiterate this argument. In QED physical ob¬ 
servables are calculated in a perturbation series of the form F (e^) = 09 - 1 - 026 ^- 1 - 
046 ^ -I- • • • where e ^ 0.3 is the charge of the electron. Now if one imagines a 
world where < 0, i.e. like charges attract, it is easy to see that the ordinary 
vacuum is unstable to the production of many electron-positron pairs which 
separate into clouds of like-charged particles. At weak coupling there is a large 
tunneling barrier to overcome because one must pay for the rest mass of the 
pairs and separate them far enough for the wrong-sign Goulomb potential to 
compensate. Thus there is a finite but exponentially suppressed rate of vacuum 
decay. A Taylor series expansion in cannot capture this non-analyticity so 
the perturbation series must be divergent. 

Similarly, the perturbation series in A diverges for the toy model considered 
here. Fade approximants are more effective because they can develop isolated 
poles in the complex A plane, however they struggle to capture the strong cou¬ 
pling behaviour at any fixed order in the approximation. Fade approximants 
are better able to capture the non-analyticities of the Borel transform, how¬ 
ever, and the widely used combination Borel-Fade approximants give a better 
global approximation. This occurs because the Fade approximated Bore trans¬ 
form has poles in the Borel plane, which lead to branch cuts when the Laplace 
transform is taken to return to physical variables. Similarly, the self-consistent 
2FI approximations develop branch point non-analyticities and approximate the 
exact Green function rather well in the entire complex A plane already at the 
leading non-trivial truncation. However, the branch cuts in the 2FI case arise 
because the 2FI Green function obeys self-consistent equations of motion, and 
is connected to the existence of unphysical solution branches. 

A question that naturally arises is: how do these methods compare? Both 
2FI and Borel-Fade methods have the ability to accurately represent non- 
analyticities of the exact theory and so out-perform other methods. However, 
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we conjecture that self-consistently derived equations of motion “know more” 
about the analytic structure of the underlying theory than do the generic Borel- 
Pade approximants. Hence we test the hypothesis that the 2PI methods should 
be more accurate that Borel-Pade and, indeed, find this to be the case, at least 
in certain regimes. 

The theory discussed here, although admittedly a toy model, also has phys¬ 
ical relevance. Independent from us, Beneke and Moch found this toy model 
as the theory governing the zero mode of scalar fields in Euclidean de Sitter 
space [18) . They performed an analysis very similar to ours, finding that a 
non-perturbative treatment is necessary and comparing 2PI and (Borel-)Pade 
resummed approximations. However, they present this analysis very briefly as 
part of a larger discussion of scalar helds in de Sitter space. Further, their com¬ 
parison of the 2PI and resummed techniques, while correct as far as we can tell, 
is not very detailed. Here we present detailed discussion of the interplay be¬ 
tween 2PI effective actions and various resummation techniques. Our use of the 
spectral function to quantify the non-analyticities present in the Green function 
in aid of this comparison is, as far as we are aware, a new aspect. 

This paper has a somewhat pedagogical flavour, and readers familiar with 
field theory, Borel summation and Pade approximants can skim through sections 
and 1^ where these topics are discussed, pausing only to pick up our notation 
and our derivation of the spectral function for the exact theory (Section 
and the Pade resummed theory (Section 3.3). In Section we compute the 
2PI effective action. Green function and corresponding spectral function for the 
theory. In Section we introduce, for the first time to our knowledge, a hybrid 
2PI-Pade scheme and show that it accelerates convergence at weak coupling 
while having the correct behaviour at strong coupling, like 2PI approximations 
but not the usual Pade method. Finally in Section]^ we draw our conclusions. 


2. Exact solution of zero dimensional QFT 

We consider a Euclidean QFT in zero dimensions, i.e. a probability theory 
for a single real variable q given by the partition function 

Z[K]=nJ dgexp (1) 

in the presence of a source K for the two point function. iV is a normalizing 
factor chosen so that Z [0] = 1. This theory has been discussed before in the 
context of exact and non-perturbative methods in field theory (see, e.g., [Hds] 
and references therein) because, despite the absence of spacetime, the theory 
possesses a perturbative expansion in terms of Feynman diagrams with the same 
combinatorial structure as more realistic theories. It was also discussed in m 
as an effective held theory for scalar held zero modes in Euclidean de Sitter 
space. 

We restrict attention to the > 0 theory since, though the m? < 0 theory 
exists and may be interesting for other purposes, it possesses no sensible weak 
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coupling limit and in zero dimensions does not give a broken symmetry phase 
anyway ^ In the absence of symmetry breaking we may omit a source term for 
q. The integral diverges for ReA < 0, but we can define the value of Z [K] at 
complex A by analytic continuation. Then Z [K\ possesses a branch cut along the 
negative A axis. Physically, this signals the instability of the negative A vacuum 
due to tunneling away from the local minimum at g ^ 0 to g ~ ±-^—6m^/A 
followed by rolling down the inverted quartic potential which is unbounded from 
below. The branch point at A = 0 means that the weak coupling perturbation 
series has zero radius of convergence. This is similar to the behaviour which 
exists in most theories of physical interest as argued by Dyson m- 

The analytic behaviour of the integrand for A > 0 is shown in Figure 
One can see that the integrand has a maximum at g = 0 and saddle points at 
g = ±im^6/X. Z [K] for the > 0 theory can be obtained from integration 
along the real axis, while integration along the imaginary axis gives the < 0 
theory. Similarly, taking the reciprocal of the integrand sends m? —>■ 

A —> — A and also reverses the colour map in Figure (changing maxima to 
minima and vice versa). From this we can see that, as expected, the partition 
function diverges for A < 0 and theory makes no sense. 

We introduce the conveniently rescaled variables k = if/m^ and p = 3m^/4A 
and obtain 


\/p(l + ^) exp (^p{l + kf^ Ki/4 (1 + , 


( 2 ) 


where Rri /4 (• • •) is a modified Bessel function of the second kind. This expres¬ 
sion is valid so long as 

R-e (v^(l-I-fc)) > 0, (3) 

which extends the definition Q to the entirety of the cut A-plane. The normal¬ 
ization factor is 

= (4) 


2 


and finally 


Z [K] = \/l + k exp (^p (1 + ~ ^ ) 


K 


1/4 (p (1 + ^) ^ 
Ki/4 (p) 


( 5 ) 


Z [K] is used to compute the expectation value of physical observables 
O (g^), by the standard trick of differentiating under the integral and then 


^Z[K] only really depends on the ratio X/m^, so the only sensible definition of weak 
coupling is A a limit which sends the entire support of the integral to g* ^ dioo if 

< 0. 

^The standard argument for spontaneous symmetry breaking, i.e. that the tunneling am¬ 
plitude between vacua tends to zero exponentially as the volume of spacetime tends to infinity, 
is clearly inapplicable in this case. 
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Figure 
A = 1. 
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and phase lb of the integrand in Q for + K = 1, 
The integrand has a maximum at g = 0 and saddles at q = ±imy/6/X. 
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removing the source: 


(O (q^)) dqO (q^) exp {-\mV - 

= Oi-2dK)Z[K]\^^,. 


We also define the generating function 


W[K]=-\nZ [K] 


( 6 ) 

( 7 ) 


and note that averages are found by taking derivatives of W. For example, 


d 


W[K] = \{q^)^^\G, 


52 

dK^ 


W[K] 


1 

4 


((?*> - Vf) 


K 


1 (^V^(4)(54_2(52^ ^ 


( 8 ) 

(9) 


where the subscript K indicates the average is taken at a fixed value of K. 
(Note that W [K] is not the connected generating function as usually defined 
because K is a, two-point source.) G and are the proper two and four point 
functions respectively. To lowest order in perturbation theory and with K = Q, 
G = llw? + 0{\) and 0^) =A-fC>(A2). 

The exact value of W [K] is easily obtained directly from the definition, 
giving 


W[K] = -^\n{l + k)- p {1 + kf-l 


— InK 


1/4 


(p(l-kfc)' 


lnKi/4 (p) . 

( 10 ) 


By direct differentiation we obtain the exact two point correlation function 


m 




K 


3/4 


(p(l + *)■ 


K 


1/4 (p (1 + ^)^) 


- 1 


( 11 ) 


or for the original (K = 0) theory, 


m~G = 4p ( 


/K3/4 (p) 


Ki/4 (P) 


-1 


( 12 ) 


Like Z^ G possesses a branch cut discontinuity from A = 0 to A = — 00 . At 
A = 0 one obtains the usual free (Gaussian) theory result G = Gq = m~‘^. In 
the strong coupling limit, A —)■ 00 , G ~ [2-d3r (|) /T (|)] A”^/^ _|_ (^ (A“^). G 
is shown in Figure from which the branch cut is obvious. This can also be 
seen in more detail in the complex A plane as shown in Figure One sees that 
not only does G possess a branch cut, it is analytic in the cut plane and is in 
fact a Herglotz-Nevanlinna function (i.e. G(A)* = G(A*) where * is complex 
conjugation). This means that G has a nice integral representation which we 
derive now in order to quantify the branch cut. 
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1.0 



Re[G] 

lm[G[A-f£]] 

lm[G[A+fe-]] 


Figure 2: G from (12) as a function of A for m = 1. The sign of the imaginary 
part depends on whether one approaches the A < 0 cut from above or below. 
Note the imaginary part is exponentially suppressed near the origin because the 
vacuum decay process is non-perturbative. 


We obtain the integral representation for G using the Cauchy formula 

where the contour G circles A in the counter-clockwise direction and avoids the 
cut. Deforming the contour to run on the circle at infinity and around the cut 
and using G (A) —)■ 0 as |A| —>■ oo, we can write the integral in terms of a spectral 
function cr (s) = [G (—s — ie) — G (—s -I- ie)] 0 (s) /2Tri = Im [G (—s — ie)] 0 (s) /tt 
where 0 (s) is the Heaviside step function, such that 


G(A) 



ds 


O'(g) 
5 -|- A 


(14) 


We find 


4V2 

= - 2 2 '^ 


i-pf 




(15) 


where {p) are modified Bessel functions of the first kindj^ a (A) is shown in 
Figure]^ 


^Note that the physical interpretation of this spectral function is unrelated to the usual 
one in field theory since, for one thing, there is no such thing as energy in zero dimensions. 
We consider a (s) a purely formal device that gives information about the analytic structure 
of G. 
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and phase 3b of G in (12) for = 


1 in the complex A 
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Figure 4: Spectral function a (A) of (151 for m = 1. 


3. Perturbation theory and re-summations 

3.1. Perturbation theory 

At small coupling one often uses perturbation theory in A which proceeds 
by expanding the exponential 


2 [K] = ^ j dgexp (m^ + K) 

" /_” i-E s {-I ("•' + 

£ s (-s=p (- 5 


n—0 

oo 


“Ssf-sM JV2“Ur(2„ + 


n—0 


n! V 4! 


(m2-hA:)^"+ = 


(16) 


In the third line we have formally interchanged the sum and integral, leading to 
an asymptotic rather than convergent series for Z [K]. Again, N is determined 
by Z [0] = 1, giving 


N = 


1 f 1 


.n—0 




1 


2',2ri+i 


(to^) 


so 


Z[K] = 


E“oi^(-irA)”2-r(2n+l) 


( m 2 + A :)^"+2 


E“oi^(-iTA)”2-r(2n+A)^^^ 


(17) 


(18) 
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and 


OO - 

HqAi = -i„^- 

n—0 

From this we find 




22nr 


2n + 


(m^ + K) 


2n+i 


f + const. (19) 


G = 2dKW 


1 

m? 

1 

w? 


4 Er=i(F^(-iTA)”22"r(2n+i);;^ 

E“o^(-iAr2-r(2n+4)^ 


A 2A2 _ 11A3 34A4 

2m® Sm^® Sm^"^ Qm^® 


( 20 ) 


The perturbative approximations Gn are the 0(A”) truncations of this series. 
The first few are shown compared to the exact G in Figurej^ Note that these are 
simply the low order Taylor series approximations to G. These approximations 
apparently converge poorly to the exact solution, and in fact we will shortly 
show that the series diverges. 

The series for G can be described in terms of Feynman diagrams by the 
following rules: 

1. Draw all connected graphs with two external lines (i.e. lines with one 
end not connected to any vertex) constructed from lines and four point 
vertices. 

2. Associate to each line a factor Go = 1/m^. 

3. Associate to each vertex a factor —A. 

4. Divide by an overall symmetry factor being the order of the symmetry 
group of the diagram under permutations of lines and vertices. 
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We illustrate these rules by giving the first few terms in G: 

--He 

+1 g 0 0 ^ i-o-S—+i 8 0 

8 8 8 



+ 


1 

12 





+ ^-+ O (A^) . 


( 21 ) 


W [K] can be written as a similar diagrammatic series in terms of connected 
vacuum diagrams (i.e., those with no external lines). 

The series (20) has the form G = (A/m^) where the coeffi¬ 

cients asymptotically obey c„+i ^ — |nc„ as n —>■ oo, thus the radius of conver¬ 
gence of the series is zero. This is consistent with the fact that we are perturbing 
around a branch point of the exact solution: no approximation of G in terms 
of analytic functions can converge at A = 0 because Z \K] is itself undefined 
for ReA < 0. The terms of the series start to increase when c„+i (A/m^) ~ c„, 
i.e. n ^ 3m^/2A, meaning the series is useful for A <?; but fails immediately 
for a moderately strong coupling A « m'^. This is typical asymptotic series 
behaviour as shown in Figure Extrapolating perturbation theory to strong 
coupling A ^ is simply impossible, although the exact solution is well be¬ 
haved there. (In fact G can be expanded as G = c„ (A/m‘^) 

displaying explicitly the branch point at A = oo.) 


5. .2. Borel summation 

We have seen that the series expansion for G diverges for all A 7 ^ 0. This is 
typical of perturbation series and usually signals some singularity of the exact 
solution for unphysical values of A. In our case, indeed, the theory does not exist 
for ReA < 0 and the exact solution possesses a branch cut on the negative A 
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Figure 5: Perturbative approximations to G up to O (A®) for m = 1, compared 
to the exact solution. 
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Figure 6: Relative error |(G„ — G) /G|of the n-th perturbative approximation 
to G for m = 1 and A = 1/4, showing the decreasing then increasing behaviour 
typical of asymptotic series. 
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axis, a feature which cannot be reproduced in any order of perturbation theory. 
However, the perturbation series is asymptotic and does contain true informa¬ 
tion about the exact solution, even if A is large enough that the series is not 
useful practically. Because of the ubiquity of this phenomenon, mathematicians 
have invented a number of series summation techniques which assign a finite 
value to certain types of divergent series and which obey certain consistency 
properties (e.g. the value assigned to a convergent series is just its naive sum). 
Here we investigate Borel summation, which is capable of summing factorially 


divergent series like (201. 


Suppose that divergent series but that the Borel transform of 

the series, defined as 


f){x) = 


OO „ 

^ anX'^‘ 


n—0 


converges for sufficiently small x. Then, if the integral 


poo 

B{x)= / e“‘(/) (xt) dt 

Jo 


( 22 ) 


(23) 


exists the Borel sum [T6l[T7l[T8] of the divergent series is defined as B — 

B{1). This definition is justified by substituting the series for (j>{xt) into the 
integral and evaluating term-wise and noting that B (x) ^ 
main drawback of Borel summation is that one must know the precise form 
of a„ for all n to compute (j){x), which is rarely the case in field theory. For 
this reason Borel summation cannot be usefully applied directly. However, one 
may use Fade approximants as discussed in the next section to recast the Borel 
transform in a useful way. 

We note that the key to Borel-summability of the perturbation series is the 
alternating sign (—1)" of the n-th order term. To see this consider the two series 


= ( 24 ) 

n—0 

OO 

,52 = ^ A”n!, (25) 

n—0 

which differ only by the alternating sign. The Borel transforms 4>i,2 (x) are 


OO ^ 

^■(-) = E< =i + A.. 

n=0 

(26) 

OO ^ 

f> 2 (x) = E(^'^) = 1_ Ax’ 

n=0 

(27) 
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and the Borel sums are 


roo t 

B[Si] = B,il) = J (28) 

e“‘ 

B[S2] = B2il) = -^dt. (29) 

In the first case the integral exists and B [S'!] = A“^e^/^r (O, where T {a, b) = 
is the incomplete gamma function. However, the second inte¬ 
gral hits a pole at t = 1/A. There is no natural prescription for avoiding the 
pole, which leads to an ambiguity in the sum of ±7r*A“^e“^/^. This is a non- 
perturbative ambiguity called a renormalon |19j . In every known case where 
this arises in field theory the renormalon is connected to a non-perturbative 
finite action solution of the field equations, i.e. an instanton or soliton, and a 
correct evaluation of the path integral which sums over all saddle points (not 
just perturbative ones) removes the ambiguity. Key to the practical application 
of Borel summation is the location and classification of all renormalons in a 
given theory m- Sophisticated techniques have been developed to deal with 
this situation which are beyond the scope of this paper [Miiniin]- 


3.3. Bade approximation 

Borel summation on its own has limited usefulness in practice because one 
often only knows a few low order terms of perturbation theory, and the potential 
existence of renormalon singularities. There exists another technique which of¬ 
ten improves perturbation series and is far more useful in practice (and is often 
combined with Borel summation). Bade approximation approximates a func¬ 
tion by rational polynomials which generally converge rapidly, are very useful for 
numerical computation and help to estimate the location of singularities of the 
function in the complex plane. Many software packages include standard rou¬ 
tines for evaluating Bade approximants, for instance the BadeApproximant 
function in Mathematica. In this section we apply Bade approximants directly 
to the Green function and then to the Borel transform. We find that the latter 
approach is clearly the better one. 

The [N, M)-Bade approximant of a function “nic” is given by [IB] 


Pm{^) 


A 

Z^n=0 

B 

Z^n=0 


(30) 


where without loss of generality one takes Bg = 1. The remaining JV + M + 1 
coefficients are chosen so that the Taylor series of (x) matches the pertur¬ 
bation series up to O . Due to the denominator, Bade approximants 

develop poles in the complex x-plane, allowing the close approximation of more 
singular functions than Taylor series are capable of. Examples of the use of Bade 
approximants in field theory can be found in [niiii] and references therein. 

Now we find the Bade approximants to m^G with x = Xjrn^. Here we 
restrict attention to the approximants where M = N + 1. This guarantees that 
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Table 1: First few Fade approximants for im?G. 


N 

r>N 


1 

0 

1 j_ 

^ ^ 2mi 

1 

1 + ^ 

1 1 5A I 7A2 

2m^ 12m8 

9 

1 1 16A 1 59A^ 

3m^ 12m« 


1 1 35A 1 43A2 I 7Va3 

' 6m‘^ 6m® ' 72m^2 

3 

1 , lOA , 155A3 , 15A3 
' ' 6m® ' m^2 

1 1 21A 1 365A2 1 295A3 , 385A<i 

2m^ 12m» ' 12mi2 ' 144mi6 

A 

1 1 16A 1 315A3 I 1190A3 , 7945A^ 

' m*^ ' 4m® ' 9 m ^2 1 


1 1 33A , 259A2 , 11935A3 , 14315A<i , 7315A5 

_ 2m.* 3m8_'_72771,13_'_144ml8_' 864m,30 


t 0 as A —> 00 . If we had used the usual diagonal approximants {N = M) 
we would find an unphysical constant term > An/Bn as A —> oo. Note 

that it is impossible to match the true 1 / a/A behaviour of G as A —> oo using 
Fade approximants centred on the origin. The best that is possible in this limit 
is ^ A“^. (As it happens, using x oc VX does not allow one to resolve this 
issue: the same approximants are found only with x'^ everywhere in place of 
X. This is because the l/VX behaviour is due to the branch point at infinity, 
which is infinitely far from the origin where the Fade approximants are matched 
to perturbation theory. Low order Fade approximants can extract information 
about the branch point near the origin, but evidently not the one at infinity.) 
The first five approximants for m^G are shown in Table and the first three 
are plotted in Figurej^with comparison to the exact G. Note that the existence 
of the integral representation (14 1 for G implies that G is a Stieltjes function, 
meaning one can prove convergence properties for the Fade approximants as 
N,M —)■ 00 , though we are not concerned with this analysis here (see [TB] for 
details). 

The {N, N + 1)-Fade approximant can also be written as 


N 




N 

N+1 


= E 

i=0 


X-Pi' 


(31) 


where and pi are the Tth residue and pole respectively. Note that since all 
of the coefficients in the denominators of Table are positive and real, all of 
the poles must either be on the negative real axis or they must be complex and 
come in complex conjugate pairs. Numerical experiments suggest that all the 
poles lie on the negative A axis, though we do not know a proof of this for all 
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Re[G] 

Pade° 

Pade2 

Padel 


Figure 7: First three Fade approximants to G. Note the simple poles developed 
for A < 0. 


N. Assuming this is generally true, Fade approximants give a representation 
of G which approximates the continuous spectral function by a sum of delta 
functions 

N 

(T (s) « ctat (s) = ^7-^(5 (s+P j) . (32) 

i=0 

As 77 —?> oo the poles become denser and fill the negative A axis, eventually 
merging into a continuous branch cut. Similarly the spectral function turns into 
a dense sum of delta functions which, when considered acting on any sufficiently 
smooth test function, smooths into a continuous function. The first few are 
shown next to the exact spectral function in Figure for comparison. 

Now we consider Fade approximation of the Borel transform of G. First we 
note the following connections between the Borel transform (j) and G and a: 


V 



— (f) (x) > s 

X^V X^S 



(33) 


That is, the Green function is related to the Laplace transform of the Borel 
transform, while the spectral function is related to the inverse Laplace transform 
of the Borel transform. These relations can be shown using the definitions of (j) 
and cr and the integral representation of the (inverse) Laplace transform. This 
allows us to extract the spectral function directly from the Borel transform. 
Note that each pole of </> yields by the inverse Laplace transform a term of the 
form A“^ exp (—fc/A) in a, where k is controlled by the location of the pole. The 
general Borel-Fade approximation for cr is a superposition of terms of this form. 

We show the first few approximants to G and a in Figures]^ and W respec¬ 
tively. One sees that the low order Borel-Fade Green functions are reasonably 
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cr(Pade°) 

cr(Pade2) 

a(Pade3) 

a(Pade4) 


Figure 8: Fade approximations to the spectral function a (A) consist of an 
increasingly dense set of delta functions. (Note the delta functions have been 
smoothed for visual purposes only.) 


accurate (within a few percent for A < 5m^) but the spectral functions are not 
approximated particularly well. Certain approximants to a oscillate erratically 
and even become negative for certain values of A. Except for the fact that the 
best approximant is the highest order one plotted, there is no clear sense in 
which the Borel-Pade approximations appear to converge to a. However, even 
this bad approximation is at least a continuous function, as opposed to a sum 
of delta functions. 

4. 2PI Approximations 

The 2PI effective action is a functional of the Green function G = ( 9 ^)^ 
which is defined by the Legendre transform 011] 

P [G] = VF [K] - KOkW [K] (34) 

where K is solved for in terms of G. P [G] obeys the equation of motion 

dGT[G] = -lK. (35) 

The standard derivation of the 2PI action 00] gives in this case 

P [G] = ^ In (G“^) + + 72 PI + const., (36) 

where — 72 P 1 is (minus due to Euclidean conventions) the sum of two particle 
irreducible vacuum graphs, i.e. those graphs which do not fall apart when any 
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Figure 9: Borel-Pade approximations to the Green function G (A) (ratio of ap- 
proximant / exact). Approximants are calculated by numerical integration of 


(23). 
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Figure 10: Borel-Pade approximations to the spectral function a (A) (ratio of 
approximant / exact). Note that the (2,2) approximant is off scale. 
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two lines are cut, where the lines are given by G and vertices by — A. Explicitly: 


- 72 PI = 



O (A^) . (37) 


The equation of motion in the absence of sources is 

G~^ =m^ + 2i9g72PI, 
which has the diagram expansion 


(38) 


G-^ =■ 


1 




1 




C>(A^ 


(39) 

Notice that there is a dramatic reduction in the number of graphs of a given or¬ 
der compared to perturbation theory ( 21 ) due to the two-particle irreducibility. 


This is one of the major benefits of 2PI approximation schemes in realistic the¬ 
ories, such as gauge theories, where the Feynman diagrams proliferate rapidly. 

The equation of motion is Dyson’s equation and the second term on the 
right hand side, —E = 29 g 72 pi, represents the exact one-particle irreducible 
self-energy of the propagator G. This can be put into the usual form of a Dyson 
equation by noting m? = Gq ^ and multiplying both sides by GGq: 


G = Go + GoEG. (40) 

Iterating this gives the infinite series 

G = Go + GoEGo + GoEGoEGo + • • • . (41) 

The main difference between this and the usual perturbative Dyson equation is 
that E contains exact propagators G rather than Gq. Inserting the expression 
above for G into E one finds that, even if one retains only a finite number of 2PI 
diagrams in 72 P 1 , E contains an infinite series of perturbative self-energy graphs. 
This is the motivation in the literature for talking about 2PI as a resummation 
method. 

By power counting (or counting line ends in the corresponding diagrams) 
we find 72 P 1 = 9 '" possible to derive the 7 „ by considering 

the symmetry factors of the two particle irreducible Feynman diagrams, but we 
can also obtain the 7 „ in a simple automated way using knowledge of the exact 
solution G = G. Substituting G into the equation of motion and the expansion 
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for 72 PI, expanding about A = 0 and matching powers of A we can determine 


_ G^A G^A^ G^A^ 5G®A4 IOIGIOA^ gSG^^A^ 

“ ^8 4^ ^8 128 960 256 

8143G14A^ 271217G1®A® 374755G1®A® 5151939G20ai° 

5376 36864 9216 20480 

697775057G22A11 3802117511G24a12 201268707239G26a13 

405504 294912 1916928 

11440081763125G28A14 , 5148422676667G30a15 1665014342007385G32ai® 

12386304 589824 18874368 

4231429245358235G34ai^ 921138067678697395G36ai® 

4456448 84934656 ^ ' 

(42) 


We do not know of any closed form expression for either the coefhcients of 
this series or its sum (implicit analytic expressions can be derived; however, 
these require the inversion of G = G (A) for A(G), which is not known to us 
in closed form). However, after the first few terms the coefficients seem to be 
well approximated by 7^+1 ~ — |z 7 i, the same as for the perturbation series. 
This has the hallmark of an asymptotic series. Like perturbation theory, the 
2PI series does not converge. 

The first non-trivial contribution to the equation of motion gives 

G(ij = +-G(i), (43) 


where the subscript (1) indicates terms of order O (A^) have been kept. This 
has two solutions 


G(i) = 


—± Vm'^ + 2A 

A ■ 


(44) 


One of these solutions is unphysical and we must choose the + sign. As A —>■ 0, 

A2 




TTi? 2 m® 2 m^® 


+ G(A3 


(45) 


which matches perturbation theory up to O (A^) terms as expected. However, 
unlike perturbation theory, the strong coupling limit A —)■ 00 exists and gives 


G(i)^ 



(2A)^/^ 


+ 0 



(46) 


This series has the correct form in powers of though the leading coefficient 

is incorrect by « 15%, and the sub-leading coefficients are incorrect by « 40%, 
70% and 100% etc. Nevertheless, it is remarkable to achieve any accuracy at 
all given the simple nature of the approximation and the fact that 72 P 1 was 
truncated at leading order in A! G(i) is a much more uniform approximation to 
G than the perturbative approximation mT^ — . 


21 




Figure 11: Comparison of exact spectral function u (A) (151 with the two loop 
2P1 approximation cr^i) (A) (48) for m = 1. Already, the simplest nontrivial 2P1 
truncation gives a much better approximation than perturbation theory or Pade 
approximants to arbitrary order. 


This result is possible because of the branch cut (7(1) possesses on the neg¬ 
ative A axis. The discontinuity across the cut 


G(i) (A + le) - G(i) (A - le) = ^^ 8 (-A - ^ 


gives the spectral function 


(7(1) (A) = 




ttA 




(47) 


(48) 


This is a far better approximation to the exact spectral function than obtained 


from any of the other techniques, as can be seen from Figure 11 


We show all of the first order approximations in Figure 12 Note that for 
A/m"^ > 0 the best approximation is the 2PI, followed by the Borel-Pade, then by 
Pade, then perturbation theory last of all. The situation for negative A is com¬ 
plicated. The best approximation overall is the 2PI, though it has an unphysical 
cusp where the branch cut starts (A/m"*^ = —1/2 in Figure [T^. It appears the 
2PI approximation trades sensitivity to the exponentially small portion of a in 
exchange for a better global approximation. The Pade approximation is good at 
small negative A but hits a pole at A/m'* = — 2 and there loses all validity. The 
Borel-Pade approximation is smooth and more accurate than the 2PI near the 
peak of G, but eventually becomes invalid, even negative, at sufficiently large 
negative A (A/m* < —5.37). 

At n-th order, (38) is a degree 2n polynomial in G which has 2n roots. 


only one of which is physical. For n = 2 there are analytical expressions for the 
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Figure 12: Comparison of (the real part of) the exact Green function G to each 
of the approximations discussed to first non-trivial order in each case. 


roots, though they are very bulky, and for n > 3 (38 1 must be solved numerically. 
Picking out the correct root for a given value of A is tricky in general and we 
leave this exercise to the reader. However, one can see on general grounds that 
truncations of (38) give a singular perturbation theory in A: at A = 0 the 
physical solution starts at Go and the spurious solutions flow in from infinity as 
inverse powers of A as A increases to finite values}^ At strong coupling the roots 
generically approach each other and one must carefully track them through the 
complex plane. We suspect that a resummation of the 2PI series would remove 
some or all of these spurious solutions, though we do not have a proof of this. 
We examine two potential methods to achieve this in the next section and find 
mixed results. 


5. Hybrid 2PI-Pade 


Since the series of 2PI diagrams is asymptotic, we may consider using a 
series summation method to improve convergence. Note that this is logically 
independent of the resummations embodied in the 2PI approximation itself. We 
may perform Pade summation of the action term 72 P 1 or the equation of motion 
term dG 72 Pi- We consider both. First, consider Pade summation of the action 
which matches the series expansion up to order N + M: 


r[G] = iln(G-i) + im 


^G- 


2^n=0 


(7) 


(AG2 




const. 


(49) 


also note that at least one of the spurious solutions (and always an odd number of 
them) are real since the coefficients in ( |38| l are real and complex solutions must occur in 
conjugate pairs. 
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The equation of motion becomes 


G-i =m^ + 2 


d 


(7) 


(AG 2 


dG 


(AG^r 
4ELo Ef=o - fc) 


?(7) 


= TO + 


En=0 ' 


(7) 


(AG 2 )" 


( 50 ) 


Multiplying by the denominator this becomes a polynomial equation of degree 
2 {N + M) as expected. This equation will have 2 {N + M) — l spurious solutions 
as does the ordinary 2PI approximation of matching order. 

To take a specific example, consider 72 P 1 up to O (A'^) and the matching 
(2,2)- and (1,3)-Pade approximants: 


72PI = 


G^A G^A^ G^A^ 5G8a4 


8 

G^A 

8 


48 

480 


48 


1 -f 


41G^A 

20 


7G^A2 

40 

G^A 


8(1 


G^A 

6 


5G^A^ 

36 


128 


+ ■ 


113G6A3\ 
432 ) 


(51) 

(52) 

(53) 


The solution resulting from any of these versions of 72 P 1 cannot be written in 
closed form, however numerical solutions are shown in Figurefor moderately 
small A (for A outside of this range, different roots become relevant). We see 
that the hybrid solutions are more accurate than the fourth order 2P1 solution, 
although the level of improvement depends strongly on the type of Pade ap- 
proximant employed. The trade-off is a loss of accuracy at large A, as shown 
in Figure for the best roots we can find. Each approximation decays with 
the correct leading A“^G power, however they differ by O (1) factors which are 
not negligible. Remarkably, the best approximation of those shown is the O (A) 
2P1 approximation. For the greater computational expense the (2, 2)-hybrid 
approximation is not noticeably an improvement at large coupling. The appar¬ 
ent clustering of the O (A'*) 2PI and (l,3)-hybrid approximations away from 
the exact value as A —>■ 00 probably stems from the minus sign of the leading 
term in i9g 72PI in these approximations, suggesting that one should not use 
approximations with this property. 

Now we try Fade summing the equation of motion: 


G-i 


= TO^ -I- 2 


(AG^) 

(AG 2 )" 


n 


(54) 


where the explicit factor of G~^ on the right hand side simply ensures the self¬ 
energy is an odd function of G as it must be. This equation of motion is of 
degree max (27V, 2M -|- 1), which for typical choices of N and M results in a 
rough halving of the number of spurious solutions. To obtain an analytic result 
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Figure 13: Comparison of exact G, first and fourth order 2PI and(2, 2)- and 
(1,3)-2PI-Pade hybrid solutions at small coupling. The (2,2)-hybrid solution 
lies almost on top of the exact solution. 
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Figure 14: Comparison of exact G, first and fourth order 2PI and(2, 2)- and 
(1,3)-2PI-Pade hybrid solutions at large coupling. 
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we consider the first non-trivial approximant, i.e. N = M = 1. The resulting 
equation of motion is 


G-i 


= m 


1 

+ - 


AG 


2 1 + iAG2’ 


(55) 


which agrees with the usual 2PI equation of motion up to terms of order O (A^), 
however it only has two spurious solutions instead of three. The physical solution 
obtained by matching at small A is 


^ If,, 36m^-\ 


(56) 


where 


X = A^ - STSA^m^ + ISAmV^ (-2A2 + lUm^ + 429\m*). 


Indeed, this matches perturbation theory up to O (A^) terms. However, the 
large A behaviour in this approximation is pathological: 


G 


hy 


2nn? 


18to^ 

A 


+ 0 


(a-3/2) 


(57) 


as A —>■ oo. This reflects the existence of an unphysical branch cut on the 
positive A axis starting at A/m^ = | (143 + 19v^) ~ 214, which Ghy has in 
addition to the expected cuts on the negative axis. The existence of this cut 
also renders the derivation of the spectral function invalid, meaning that Ghy 
cannot be written in the form (14|. It is not clear at this stage whether other 
forms of Fade summed equations of motion lack these pathologies. We leave 
further investigation of hybrid approximation schemes to future work. 


6. Discussion 

Two-particle irreducible effective actions are the subject of a rich litera¬ 
ture and have found applications in diverse areas from early universe cosmol¬ 
ogy to nano-electronics (e.g. [531 [33] ). The virtues of approximation schemes 
built on 2PI effective actions are often explained in terms of a resummation 
of an infinite series of perturbative Feynman diagrams which are encapsulated 
in the non-perturbative Green function G, from which the 2PI diagram series 
is built. However, the 2PI effective action is not a resummation scheme: it is 
a self-consistent variational principle. The definition of the 2PI effective ac¬ 
tion in terms of the Legendre transform is crucial for the self-consistency of the 
scheme. The immediate practical consequence is that any modification of the 
2PI effective action which does not derive from a consistent modified variational 
principle is very likely to be inconsistent. So for instance, the consistency of re¬ 
cent attempts to improve the symmetry properties of 2PI effective actions [3S| is 
guaranteed by the existence of a suitable constrained variational principle, how¬ 
ever, ad hoc attempts to modify the equations of motion to satisfy symmetries 
will fail. 
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In this work we have pointed out the distinction between resummation and 
self-consistent approximations using an exactly solvable zero dimensional “field” 
theory. The perturbation theory has zero radius of convergence due to a branch 
cut on the negative coupling (A) axis, a fact which is invisible to perturbation 
theory. The theory is Borel summable, with the Borel sum giving the exact 
answer. However, Borel summation alone is not usually this helpful in prac¬ 
tice. Bade approximants well describe the Green function at weak coupling, 
though not at strong coupling. However, the combination of Borel and Bade 
approximation yields an effective global approximation scheme for the Green 
function. 

The two point Green function of this theory admits a nice integral represen¬ 
tation in terms of the spectral function (i.e. discontinuity of the branch cut) 
which allows us to see that the Bade approximants improve perturbation theory 
by allowing the spectral function to be approximated as a sum of delta functions 
and the Borel-Bade method gives a continuous, albeit erratic and inaccurate, 
approximation to the spectral function. 

The 2BI approximation scheme surpasses perturbation theory. Bade and 
Borel-Bade approximants already at the leading non-trivial truncation. Like 
the Borel-Bade method, 2BI approximations can develop branch points and 
represent the spectral function by a continuous distribution. However, the 2BI 
approximation is quantitatively superior at the leading truncation. 

We speculate that this is because the Borel-Bade method is a widely appli¬ 
cable general “black box” method, however the self-consistent 2BI equations of 
motion are derived within a particular field theory of interest. This gives the 
2BI method “insider information,” from which it should be able to construct a 
better approximation. This comes at the cost of spurious solutions which must 
be eliminated and the added difhculty of finding the 2BI effective action in the 
first place. 

Finally we introduced (for the first time, to our knowledge) a hybrid 2BI- 
Bade scheme, using Bade approximants to partially resum the 2BI diagram 
series. The quality of the result depends strongly on the type of Bade ap- 
proximant used, with the best result in our case for the diagonal approximant. 
This hybrid approximation performs considerably better than the comparable 
2BI approximation at weak coupling, though not noticeably better at strong 
coupling. 

During proofing we became aware of a work by Kleinert |26j in some ways 
similar to our own. In [26j , Kleinert generalizes and reformulates the Feynman- 
Kleinert variational perturbation theory and compares it to the 2FI effective 
action (called by the non-standard name “bilocal Legendre transform”) for the 
toy model discussed here. He finds that the variational perturbation theory out¬ 
performs the 2FI method and notes especially the failure of the 2FI series to 
converge uniformly for all couplings (this is not a surprise in light of the fact that 
the 2FI series is asymptotic). Our work differs in scope from his in two major 
ways. First, we compare the 2FI method to traditional resummation methods 
and not to other variational methods. Second, Kleinert focuses on the value 
of the effective action itself evaluated at its extremum, while we focus on the 
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correlation function. In particular, we have focused on the analytic structure 
of the correlation function in the complex coupling plane, and introduced the 
spectral representation to aid this discussion. A study of the behaviour of 
these quantities in variational perturbation theory is an intriguing prospect for 
future work. (An intermediate result of Kleinert’s work is directly applicable to 
our discussion around (42): equation 110 of |26j gives a non-linear recurrence 
relation for the coefficients of the 2PI self-energy.) 

There are several other natural directions for extension of this work. One 
would be the calculation and comparison of higher orders in the 2PI and (Borel- 
)Pade approximations, although it is doubtful what new insights might come 
from this. It would be straightforward to extend this work to consider 4PI effec¬ 
tive actions, which depend on a self-consistent vertex function 141"^) in addition 
to gQ However, the extension to higher order nPI effective actions requires 
the introduction of new terms in the exponent of ([^ which makes the prob¬ 
lem no longer exactly solvable. A natural direction to pursue is Borel-Pade 
summation of the 2PI generating functional itself. This opens the possibility of 
eliminating spurious solutions to the 2PI equations of motion and enhancing the 
sensitivity of the 2PI method to the exponentially small regions of the spectral 
function, hence eliminating unphysical cusps. We are planning an investigation 
of these themes for the quantum anharmonic oscillator as a stepping stone to 
more physically interesting field theories. 
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